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ABSTRACT 


We  define  the  curvature  spectrum  as  the  two-dimen¬ 
sional  Laplacian  of  a  frequency-wavenumber  spectrum,  the 
derivatives  being  taken  in  the  wavenumber  directions. 
Results  show  that  the  method  has  the  characteristics  of 
a  high-resolution  f-k  estimator,  i.e.,  the  spectral  peaks 
are  sharper  than  in  ordinary  f-k  spectra.  The  curvature 
of  an  array  response  function  can  give  additional  infor¬ 
mation  about  the  array  capability,  in  addition  to  the 
array  response  function  itself.  In  particular,  defining 
the  differential  array  response  as  the  curvature  of  the 
array  response  at  k  =  0 ,  we  show  that  this  differential 
response  predicts  the  angular  resolution  capability  of 
an  array.  We  conclude  that  if  two  arra>s  have  differ¬ 
ential  responses  and  D2,  then  their  angular  resolu- 

1  '2 

tion  capabilities  0^  and  02  are  related  by  D^'  0^  = 

D3/2  0 2 .  It  is  shown  that  using  the  N-4  correlator  signal 
analysis  algorithm,  the  azimuth  resolution  of  four-  and 
five-element  arrays  with  apertures  as  small  as  20  km  is 
about  one-tenth  of  a  degree.  The  capabilities  of  some 
microbarograph  arrays  in  current  use  are  compared. 
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INTRODUCTION 


An  important  consideration  in  the  design  of  arrays 
is  the  resolution  capability  for  determining  signal 
velocity  and  arrival  azimuth,  given  a  particular  config¬ 
uration  of  elements.  This  is  conventionally  defined  as 
the  width  of  the  central  peak  of  the  array  response  func¬ 
tion  in  wavenumber  space  at  a  specified  number  of  decibels 
down  from  the  peak  gain.  For  example,  the  width  of  the 
central  peak  of  the  13-element  LAMA  microbarograph  array 
(Figure  1)  at  -3  db  is  approximately  0.02  units  of  wave- 
number.  In  this  report,  we  will  always  choose  k  to  be 
inverse  wavelength,  so  its  units  are  km"1. 

It  is  apparent  that  such  a  specification  is  impre¬ 
cise,  since  it  depends  on  an  arbitrary  gain  level,  and 
when  used  for  detection  analysis,  also  on  the  sensitivity 
and  resolution  of  the  detection  algorithm.  Furthermore, 
it  tells  nothing  about  the  shape  of  the  central  peak 
above  the  chosen  gain  level.  It  is  conceivable  that  two 
different  array  configurations  could  give  array  response 
functions  having  the  same  width  of  the  central  peak  at 
a  certain  gain  level  but  having  radically  different 
shapes,  as  illustrated  schematically  in  Figure  3. 

In  this  study  we  propose  another  measure  of  array 
resolution  capability  differential  quantity  that  depends 


only  on  the  shape  of  the  central  peak  at  its  maximum 
height,  and  compare  this  with  other  standard  measures. 
We  also  present  data  indicating  that  the  detection 
sensitivity  of  the  N4  multichannel  correlator  algo¬ 
rithm  (Brown,  1962;  Smart  and  Flinn,  1971)  is  of  the 
order  of  millibels,  not  decibels,  for  the  four-  and 
five-element  ESSA  microbarograph  arrays. 


THEORETICAL  DISCUSSION 

The  usual  expression  for  frequency-wavenumber  spectra 
(of  which  the  array  response  function  is  a  special  case)  is: 

N  N 

P(w,k)  =  7  y  S . . (w)exp[ik* (x .-x. ) ]  (1) 

1=1  3=1  J 

Our  method  is  to  compute  the  two-dimensional  Laplacian: 

vk  =  V7k  <2) 

of  the  f-k  spectrum  in  wavenumber  space: 

o  N  N  ? 

VjP(fc>,k)  =  -  l  l  Si.(a))|xi-x,  |  exp[ik*(x.-x.  )]  (3) 

i=l  ]  =  1  J  J  J 

This  quantity  is  also  a  function  of  w  and  k,  which  is 

related  to  the  shape  of  curvature  of  the  f-k  spectrum 

surface . 

When  applied  to  the  array  response  function 
'  N  N 

A(k)  =  l  l  exptik* (x.-x. )]  (h) 

i=l  j=l  ~  -1 


and  evaluated  at  the  origin  (i.e„,  the  point  k  -  0)  ,  this 
is : 


D(k)=7£A(k)  = 


N  N 

-  I  l  IXi-xJ 

i=l  j=l 


(5) 
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We  define  this  as  the  differential  array  response ,  since 
it  is  not  a  function  of  k,  but  depends  only  upon  the 
coordinates  of  the  array  elements. 

For  cylindrically  symmetric  surfaces,  i.e.,  surfaces 
which  are  symmetric  about  the  frequency  axis  k  =  0,  the 
two-dimensional  Laplacian  at  the  peak  is  also  related  to 
the  radius  of  curvature  of  a  sphere  passing  through  the 
surface  at  that  point.  It  is  because  of  this  relation 
that  we  use  the  term  "curvature  spectra". 
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RESULTS 


The  results  of  this  study  are  given  in  Table  III  and 
the  figures.  As  an  example  of  real  data,  we  chose  the 
acoustic-gravity  wave  signals  generated  by  an  atmospheric 
nuclear  explosion,  recorded  by  the  array  of  ESSA  micro¬ 
barographs  at  LAMA.  Information  about  the  signal  is  given 
in  Table  I;  Table  II  gives  the  coordinates  of  the  micro- 
barograpii  array  elements.  For  a  description  of  instrumenta¬ 
tion,  see  Cook  and  Bedard  (1971). 

These  contour  plots  presented  here  were  computed  by 
the  method  described  by  McCowan  (1968) .  The  functions 
plotted  are  described  in  the  figure  captions.  For 
simplicity  of  comparison  with  the  other  spectra,  the 
negative  of  the  curvature  spectra  are  plotted. 

Figure  4  shows  the  frequency-wavenumber  spectra  of 
the  acoustic-gravity  wave  signal  at  four  periods; 

Figure  5  shows  high-resolution  f-k  spectra  (McCowan  and 
Lintz,  1968)  at  the  same  four  periods;  and  Figure  6 
shows  curvature  spectra,  again  at  the  same  four  periods. 

Figure  2  shows  the  curvature  array  response  of 
the  thirteen  innermost  elements  of  LAMA,  over  the  same 
wavenumber  range  as  the  ordinary  array  response  in 
Figure  1.  It  is  apparent  that  the  central  peak  is 
narrower  by  about  half  in  the  curvature  response.  Thus 
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we  should  expect  a  "high-resolution"  effect  (McCowan 
and  Lintz,  1968)  in  the  curvature  spectra.  In  addition, 
the  curvature  response  enables  us  to  see  more  detail 
in  the  regular  array  response.  For  example,  the  ears 
of  the  -9  db  contour  in  the  regular  response  appear  as 
small  side  lobes  in  the  curvature  response. 

The  high- resolution  spectra  were  computed  with  a 
signal- to-noise  ratio  parameter  0.01  (McCowan  and  Lintz, 
1968)  . 

Comparing  these  results,  we  see  that  the  curvature 
spectra  show  a  resolution  greater  than  the  ordinary  but 
not  as  high  as  the  high-resolution  f-k  spectra.  It 
appears  that  at  128  seconds  period  the  high-resolution 
method  resolves  two  peaks,  where  neither  of  the  other 
two  methods  resolves  more  than  one  peak. 

Figures  7,  8,  and  9  show  the  ordinary  Huancayo  array 
response  contoured  in  decibels,  in  centibels  (equal  0.1 
to  decibel),  and  millibels  (equal  to  0.01  decibel).  It 
is  significant  to  note  that  as  the  contour  levels  become 
closer  and  closer  to  the  peak  value  (Figure  9),  the 
central  peak  becomes  nearly  circular .  Thus  it  appears 
that  we  are  justified  in  representing  the  peak  of  the 
central  lobe  of  an  array  response  as  a  spherical  surface 
in  the  region  around  the  peak  value. 


Array  resolution  capability 

Figure  10  shows  a  segment  of  output  from  our  N-4 
correlator  program  using  as  input  an  artificial  set  of 
records  containing  a  monochromatic  plane  wave  with  period 
20  seconds,  velocity  300  m/sec,  and  arrival  azimuth  180°, 

recorded  at  the  innermost  five  elements  of  the  LAMA  arrav 

✓  3 

Four  beams  were  used,  each  with  velocity  300  m/scc,  and 
with  azimuths  spaced  one-half  degree  away  from  the  known 
arrival  azimuth .  It  is  obvious  that  this  array,  whose 
effective  aperture  would  be  about  30  km,  judging  from 
the  half-power  points  of  the  array  response,  is  theo¬ 
retically  capable  of  resolving  azimuth  shifts  of  about 
one-tenth  of  a  degree,  to  take  a  conservative  lower  limit. 
This  result  is  consistent  with  results  obtained  using 
the  analog  version  of  the  N-4  correlator  on  four-  and 
five- element  HSSA  microbarograph  arrays.  We  note,  how¬ 
ever,  that  this  resolution  is  attained  on  the  analog 
version  of  the  correlator  only  by  "hand-setting"  the 
device,  since  in  the  normal  continuous  beamsweeping 
generating  mode  there  is  a  smearing  of  azimuth  due  to 
the  fact  that  the  beam  direction  moves  while  a  given 
segment  of  data  is  in  the  tape  loop  (11.  Mat  he  son ,  personal 
communication,  1909).  The  digital  version  of  the  corre¬ 
lator  does  not  suffer  from  this  drawback,  but  clearly 
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because  of  the  finite  sarupling  interval  there  is  a  mini¬ 
mum  azimuth  resolution  for  a  given  size  of  the  array  (see 
below) „ 

The  period  and  phase  velocity  20  seconds  and  300 
in/sec  correspond  to  a  region  of  f-k  space  centered  at 
wavenumber 


J_  _  1 

cT  ~  6 


0.167  km 


(10) 


The  width  of  a  region  in  wavenumber  space  centered  at 
this  wavenumber  and  intercepting  one  degree  of  azimuth 
is : 


k  =k0 


1 

cT 


7T 

180 


2 . 9 • 10" 3 


(11) 


Referring  to  Figure  9,  we  can  see  that  this  represents  a 
difference  in  gain  of  less  than  one  millibel  at  the 
lluancayo  array „  A  more  precise  measurement  gives  the 
results  0.88  millibel  in  the  east-west  direction  and 
0.99  millibel  in  the  north-south  direction.  Thus  we  arc 
able  to  say  that  a  human  operator  working  with  the 
analog  correlator  under  these  conditions  can  probably 
distinguish  gains  which  differ  by  about  0.95  millibel. 
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We  pointed  out  that  the  two-dimensional  Laplacian 
of  a  spherical  surface  is  inversely  proportional  to  its 
radius  of  curvature .  For  a  spherical  surface  the  radius 
of  curvature  is  proportional  to  the  width,  so  we  can 
say  that  if  two  arrays  have  differential  responses  I), 
and  T>2 ,  then  their  angularular  resolution  is  related  by: 


(12) 


Using  equation  (12)  and  knowing  the  angular  resolution  of 
an  array  which  uses  some  given  detection  algorithm,  we 
can  calculate  the  angular  resolution  of  any  other  array 
under  the  same  conditions  using  the  same  detection  algo¬ 
rithm,  This  is  done  in  Table  III,  taking  0„1  degree  for 
the  five-element  LAMA  array  as  a  standard.  An  analog 
correlator  processing  the  thirteen-element  LAMA  array 
could  theoretically  resolve  arrival  azimuths  to  within 
about  one  minute  of  arc. 

I  his  analysis  is  true  only  for  the  central  peak  of 
the  array  response.  If  the  array  is  mis-steered  badly 
enough  so  that  a  sidelobe  peak  is  being  seen,  then  the 
results  are  obviously  inapplicable.  This  is  an  unfortunate 
consequence  of  a  differential  quantity:  it  tells  nothing 
about  the  other  peaks.  However,  ure  will  assume,  as  those 
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who  operate  the  analog  correlator  do,  that  we  know  in 
advance  where  in  wavenumber  space  to  expect  infrasonic 
signals,  i0c0,  at  sonic  velocities  and  wavenumber.  In 
practice  this  difficulty  arises  only  at  the  shorter 
periods.  This  theoretical  capability  is  also  obviously 
degraded  if  the  signal  waveform  changes  across  the  width 
of  the  array,  or  if  there  is  significant  noise  in  the 
data, 

hffcct  of  sampling  rate  on  beaming  accuracy 

A  finite  sampling  rate  imposes  a  limit  on  azimuthal 
resolution :  for  a  given  arrival  direction,  we  obviously 
cannot  tell  the  difference  between  two  beams  so  close 
together  that  the  two  sets  of  beam  shift  indices  arc 
the  same  integers.  To  investigate  whether  this  is  a 
limitation  on  the  digital  N-4  correlator  processor  opera¬ 
ting  on  arrays  of  the  size  of  the  present  microbarograph 
arrays,  we  used  our  artificial  monochromatic  plane  wave 
data,  and  searched  over  an  azimuth  grid  so  fine  that 
the  beam  shift  integers  differed  from  one  another  by 
only  one  integer  at  one  of  the  five  channels.  The  azimuth 
intervals  were  0,5° „  For  example,  the  beam  toward  17S„0° 
had  shift  integers  45,  0,  -19,  59,  -40;  the  beam  toward 
178.5°  had  shift  integers  45,  0,  -20,  59,  -40„  Figure  19 


"r-  1 WKPq 


shows  a  comparison  of  two  outputs  of  these  runs,  one  for 
a  sampling  rate  of  4  Hz  and  one  for  a  sampling  rate  of 
1000  Hz,,  It  is  clear  that  at  sampling  rates  at  least  as 
low  as  4  llz  the  finite  sampling  interval  does  not  degrade 
the  correlator  output. 
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CONCLUSIONS 


Curvature  spectra  convey  roughly  the  same  information 
as  ordinary  or  nigh-resolution  f-k  spectra,  with  resolu¬ 
tion  intermediate  between  the  two. 

The  curvature  of  an  array  response  function  contains 
inf 01 mation  about  an  array  not  supplied  in  the  array 
response  function  itself.  For  example,  the  polarity  of 
side  lobes,  i.e.,  whether  they  are  small  peaks  or  valleys, 
is  better  judged  from  curvature  response  plots. 

1  he  differential  response  of  an  array  provides  a 
means  of  comparing  the  angular  resolution  of  different 
arrays  having  the  same  signal  detection  sensitivity. 
Azimuth  resolution  of  the  infrasonic  analog  correlator 
detector  is  0.95  mb.  At  20  sec.  period  and  200  m/sec 
phase  velocity,  experience  shows  that  the  angular  resolu¬ 
tion  of  the  Huancayo  microbarograph  array  is  about  one 
degree  of  azimuth,  which  implies  that  theoretically  the 
angular  resolution  of  LAMA,  measured  by  the  same  tech¬ 
nique  is  about  one  minute  of  arc. 
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1 

tabu; 

1* 

l.vent  Data  for  the  Acoustic 

-Gravity  Wave  Signal 

j 

Array 

LAMA 

j 

Number  of  elements 

]  3 

Sampling  rate 

2.0  samples /sec 

j 

Time  window 

02:25:44  -  20:34:16  GMT 
(25  August  1968) 

Group  velocity  window 

394  -  384  m/sec 

Date 

24  August  1968 

Latitude 

22S 

1 

Longitude 

1  391V 

| 

Seismic  magnitude  (m,  ) 

5.0 

i 

Lpicentrai  distance 

8300  km 

Azimuth  of  epicenter 
from  array 

211°  east  of  north 

*  Based  on  information  taken  from  the  NOAA/NOS  Preliminary 
Determination  of  Epicenters  releases. 


TABLE  1 1 


Array  Element  Coordinates 


Element  X 

B1 
AO 
B3 
C4 
B4 
Cl 
C2 
B2 
C3 
D3 
D4 
D1 
D2 

Element 

C7 
C9 
Cll 
CJ  3 


LAMA 

Coordinate  (km) 
9 .98 
0.00 
-7.21 
-11.59 
-1.59 
7.  23 
16.05 
4.55 
-2.11 
-19.79 
-12.27 
23.39 
16.23 

Hvancayo 
Coordinate  (km) 

0 .00 

2.4  8 

7.4  3 
4  .  36 


y  Coordinate  (km) 
7.08 
0.00 
-3, 27 
5.  20 
8.82 
16.77 
-2.11 
-5.96 
-12.72 
-15.41 
28.12 
16.85 
-20.61 

y  Coordinate  (km) 
0.00 
-4.12 
-2.08 
3.65 


x 


- . - . . 


Table  III 

Differential  Responses  for  Several  Microbarograph  Arrays 


Array 

Number  of 
Elements 

Differential 

Response 

Predicted  Angular 
Resolut  icr.J’: 

(square  km) 

(degrees ) 

Penas 

4 

_4  o'? 

0.26 

Huancayo 

4 

-500 

0.26 

LAMA 

C 

'inner  thirteen 

elements)  13 

-1.11x10 

0.018 

Washington 

5 

-1520 

0.15 

Washington  (without  DALE)  4 

-3  050 

0.18 

LAMA  (inner  five 

elements )  5 

-3402 

0.1  (assu:. 

^Assuming  a  300  m/sec  wave  at  20  sec  period. 
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Figure  1.  LAMA  array  response;  contours 
shown  for  128  seconds  period. 


Figure  10.  N4  correlator  output  for  an  input  consisting  of  a  mono¬ 
chromatic  plane  wave,  20  seconds  period,  velocity  300  meter/second 
(a)  sampling  rate  infinite;  (b)  sampling  rate  4  Hz. 


